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Summary 

Turbine  cooling  system  engineers  often  evaluate  designs  using  experiments  that  measure  heat 
transfer  rates  and  pressure  drop.  Engine  Reynolds  numbers  are  reproduced  in  large-scale  test 
sections  and  the  dimensionless  heat  transfer  coefficients  then  relate  directly  to  engine  values. 
The  experimental  test  strategy  allows  the  designer  to  quantify  the  performance  of  new,  for 
example  more  intricate,  flow  schemes.  One  common  method  entails  the  use  of  a  model 
manufactured  from  a  thermally  insulating  material  which  is  subjected  to  a  rapid  change  in 
inlet  gas  temperature.  This  paper  is  concerned  with  a  method  of  processing  these  tests  by  only 
analysing  the  gas  temperature  changes.  The  approach  offers  data  at  lower  resolution 
compared  to  conventional  methods  but  can  provide  information  in  situations  where  full 
surface  temperature  changes  can  not  be  measured.  The  advantages  and  disadvantages  of  the 
method  are  discussed  and  results  are  compared  to  the  data  from  a  conventional  analysis  of 
liquid  crystal  coated  models. 

Background  and  purpose  of  research 

Several  experimental  methods  have  been  used  in  the  evaluation  of  turbine  cooling  system  heat 
transfer  performance  by  numerous  researchers.  The  techniques  used  include  the  steady  state 
conducting  element  method  (Florschuetz  et  al.,  1981)  and  the  steady  state  uniform  heat  flux 
approaches  (Hippensteele  and  Russell,  1984).  A  transient  heat  transfer  method  using  a 
temperature  sensitive  paint  was  first  applied  to  a  turbine  internal  cooling  system  by  Clifford  et 
al.  (1983).  Since  this  early  paper,  there  have  been  numerous  examples  of  the  technique  to 
blade  cooling  research.  Baughn  et  al.  (1989)  compared  the  transient  to  the  uniform  heat  flux 
method.  Recent  transient  method  publications  include  Son  et  al  (2000),  Azad  et  al.  (2000) 
and  Gillespie  at  al.  (2000).  For  a  detailed  description  of  the  method,  the  reader  is  referred  to  a 
review  of  the  technique  given  in  Ireland  and  Jones  (2000).  The  technique  benefits  from 
reasonably  straightforward  model  manufacture  and  is  able  to  provide  high-resolution  data 
from  short  duration  experiments.  Figure  1  shows  a  rig  used  at  Oxford  to  measure  rib  cooling 
system  performance.  One  advantage  of  the  standard  transient  method  is  that  local  heat 
transfer  coefficients  are  determined  over  any  region  of  interest  for  which  surface  temperature 
data  is  available. 

The  analysis  route  covered  in  this  paper  offers  data  at  much  lower  resolution,  but,  under  some 
circumstances,  the  area-averaged  data  can  be  of  great  use.  The  reasons  behind  the 
development  of  the  technique  are  listed  below. 

1.  Independence  from  surface  temperature  thermometry.  The  technique  does  not 
require  the  surface  temperature  to  be  logged  during  the  experiment.  The  result  is  that 
experiments  can  be  performed  for  rigs  for  which  optical  access  is  not  possible.  This 
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can  be  useful  in  complex  cooling  systems  for  which  it  is  impossible  to  gain  optical 
access  to  the  full  heat  transfer  surfaces  because  of  a  limited  number  of  cameras  or 
potential  sites  for  cameras.  Often,  the  lighting  requirements  for  one  surface  prevent 
another  surface  from  being  properly  illuminated  or  viewed.  The  finite  thickness  of  an 
encapsulated  liquid  crystal  coating  limits  the  maximum  heat  flux  that  can  be 
measured  in  practice  -  see  Ireland  and  Jones  (2000).  The  new  technique  opens  the 
way  for  higher  heat  flux  experiments  to  be  performed. 

2.  Reduced  manufacture  time.  The  tests  analysed  using  the  new  procedure  have  used 
existing  Perspex  models  designed  for  conventional  transient,  heat  transfer 
experiments  using  liquid  crystals.  However,  analysis  of  only  gas  temperature 
changes  means  that  model  requirements  for  optical  access  can  be  relaxed.  This 
means  that  the  cooling  system  model  can  be  produced  from  material  that  allows  more 
rapid  manufacture  than  the  Perspex.  A  softer  plastic  material  such  as  a  wax  is  easier 
to  machine  to  the  required  scale  size.  The  material  does  not  need  to  be  transparent 
and  does  not  need  to  be  polished  as  the  heat  transfer  experiments  rely  on  gas 
temperature  measurements.  The  model  is  fitted  with  gas  thermocouples,  as  before, 
though  black  ink  and  liquid  crystal  is  not  required. 

3.  Reduced  rigging  time.  The  model  can  be  fitted  to  the  same  rig  used  for  transient  heat 
transfer  testing.  However,  there  is  no  need  to  position  lights  or  to  align  cameras  to 
illuminate  and  observe  the  liquid  crystal.  The  rig  set  up  time  is  considerably  reduced 
and  the  need  for  a  crystal  calibration  stage  is  also  eliminated. 


Figure  1  Apparatus  used  for  rib  heat  transfer  experiments. 

Application  of  gas  temperature  measurements  in  transient  experiments 

The  heat  transfer  coefficient  distribution  in  a  conventional  transient  heat  transfer  experiment 
is  calculated  from  a  procedure  that  employs  a  functional  relationship  between  the  local 
surface  temperature  change  and  the  local  heat  transfer  coefficient.  Knowledge  of  the  relevant 
gas  temperature  is  required  to  determine  the  h,  together  with  the  model  initial  temperature, 
crystal  calibration  and  thermal  properties  of  the  model  material.  The  method  requires  the  heat 
transfer  driving  gas  temperature,  specifically  Tg  in  the  equation  (1)  for  local  heat  flux,  q, 

q  =  h(rg-Ts )  (1) 

to  be  measured  directly  or  derived  from  experimental  data.  A  decision  has  to  be  made  on  what 
driving  temperature  to  use  for  the  calculation  of  heat  transfer  coefficient.  In  the  simplest  case, 
the  hs  are  calculated  based  on  an  upstream  driving  temperature  measured  using  one  or  more 
gas  thermocouples.  This  strategy  was  used  by  Ireland  and  Jones  (1985)  to  measure  the  heat 
transfer  coefficient  round  an  isolated  pedestal.  The  local  htc  maps  were  based  on  the 
centreline  temperature  of  the  developed  channel  flow  ahead  of  the  pedestal.  In  longer 
passages,  it  is  advantageous  to  calculate  the  hs  based  on  a  local  fluid  temperature.  Clifford  et 
al  (1980)  measured  the  duct  centreline  temperature  through  models  of  complex  turbine 
cooling  passages.  In  cooling  systems  with  side  bleed,  eg  Gillespie  et  al  (2000)  and  Hwang  et 
al  (1999),  it  can  be  difficult  to  establish  a  representative  heat  transfer  driving  temperature. 
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After  the  htc  data  have  been  acquired,  in  subsequent  component  analysis,  changes  in  mixed 
bulk  temperature  are  directly  related  to  total  heat  flows  so  hs  based  on  local  mixed  bulk 
temperatures  can  be  especially  useful.  There  are  also  good  physical  reasons  why  h  based  on  a 
local  temperature,  such  as  local  mixed  bulk,  are  less  sensitive  to  thermal  boundary  conditions. 
A  method  of  calculating  hs  based  on  mixed  bulk  temperature  was  introduced  by  Metzger  and 
Larson  (1986)  for  flow  round  a  bend.  A  review  of  different  strategics  was  included  in  Chyu  et 
al.  (1997)  and  later  in  Tsang  et  al.  (2000).  Gillespie  et  al.  (1994)  were  able  to  use  knowledge 
of  the  heat  transfer  coefficient  based  on  the  local  mixed  bulk  temperature  at  the  exit  to  a 
branched  duct  to  process  liquid  crystal  heat  transfer  experiments  for  which  the  entry  mixed 
bulk  temperature  was  unknown. 

In  the  present  paper,  the  surface  temperature  is  not  used  to  determine  h.  Instead,  the  gas 
temperature  signals  are  processed  in  one  of  two  ways.  The  choice  depends  on  practical 
concerns  such  as  how  rapidly  the  air  supply  temperature  can  be  switched  and  the  time 
available  for  analysis.  Both  methods  are  described  and  compared  below. 

A.  Initial  temperature  distribution  method 

For  experimental  rigs  that  can  achieve  an  abrupt  change  in  the  supply  gas  temperature,  the 
turbine  cooling  system  is  analysed  by  extrapolating  the  measured  gas  temperatures  back  to  the 
beginning  of  the  experiment.  At  the  start  of  the  experiment,  the  model  temperature  is  uniform 
and  known.  This  allows  an  energy  balance  to  be  performed  to  determine  the  intervening, 
mean  heat  transfer  coefficient. 

The  change  in  the  gas  temperature  between  planes  is  analysed  to  yield  the  number  of  transfer 
units  (defined  below)  between  the  measurement  planes.  The  NTU  is  based  on  the  average 
heat  transfer  coefficient  between  the  planes.  In  the  absence  of  bleed  flow  to  films,  the 
mathematics  is  straightforward.  For  an  initial  wall  temperature  Tw  and  initial,  mixed  bulk  gas 
temperature  changes  between  two  planes  from  Tmb(i-i)  to  Tmbi,  then 

vrt/=7A  =  ln 

me  p 

where  h  is  the  average  heat  transfer  coefficient  between  planes  i- 1  and  i,  A  is  the  wetted  area 
between  the  planes  and  m  and  cp  are  the  mass  flow  rate  and  specific  heat  respectively. 

The  initial  gas  temperature  distribution  needs  to  be  accurately  determined  for  the  initial 
temperature  distribution  method.  Attention  needs  to  be  paid  to  the  model  material  thermal 
product1  and  to  the  gas  temperature  instrumentation  bandwidth.  One  approach  is  to  ensure 
that  the  internal  surface  temperature  does  not  change  significantly  within  the  longer  of  (a)  the 
flow  establishment  time  or  (b)  the  temperature  transducer  time  constant.  Figure  3  shows  how 
long  it  takes  the  substrate  temperature  difference  to  reach  5%  of  its  asymptotic  value  under 
circumstances  in  which  the  gas  temperature  has  undergone  a  step  change.  The  authors  have 
used  Perspex  models  that  change  temperature  relatively  rapidly.  They  found  it  necessary  to 
extrapolate  the  gas  temperature  signals  back  to  the  start  of  the  experiment. 

The  tests  discussed  below  used  a  machined  finish  for  which  the  surface  average  roughness 
was  judged  to  correspond  to  a  hydraulically  smooth  condition.  It  is  worth  noting,  however, 
that,  since  there  is  no  need  to  apply  a  liquid  crystal  coating,  it  is  possible  to  machine  the 
surface  with  a  roughness  at  scale  that  is  representative  of  a  cast  blade. 

B.  Gas  temperature  variation  with  time  method 

The  gas  temperature  data  throughout  the  test  can  be  compared  to  gas  temperatures  predicted 
from  an  analytical  model  of  the  experimental  gas  temperature  rise.  The  latter  analysis 
requires  knowledge  of  the  average  h  between  planes  to  predict  the  temperature  field  so  a 
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regression  procedure  is  used  to  determine  the  most  likely  h.  The  advantage  of  this  approach 
is  that  it  can  be  used  for  experiments  where  it  is  not  possible  to  produce  a  sufficiently  abrupt 
start.  The  analysis  method  is  described  in  Tsang  et  al.  (2000)  and  is  related  to  the  approach 
introduced  by  Wolfersdorf  et  al.  (1997)  who  introduced  an  analytic  model  that  enabled  mixed 
bulk  heat  transfer  coefficient  to  be  determined  from  surface  temperature.  Similar  strategies 
have  been  used  to  assess  the  performance  of  heat  exchanger  matrices-  see  Smith  (1996). 

In  the  work  used  in  the  evaluation  of  the  new  methods,  the  temperature  at  the  exit  to  the 
ribbed  section  was  compared  to  that  predicted  using  a  numerical  model  of  the  duct.  The 
measured  inlet  temperature  was  used  as  an  input,  together  with  the  stalling  model  temperature 
and  the  measured  mass  flow  rate.  The  duct  is  dividing  into  ten  sections  and  the  wall 
temperature  and  fluid  temperature  variations  are  calculated  in  an  iterative  procedure  for  a  time 
matching  the  duration  of  the  experiment.  Regression  was  then  used  to  determine  the  heat 
transfer  coefficient  that  gave  the  best  fit  to  the  measured  exit  temperature. 


Requirement 

A.  Initial  temperature 
distribution  method 

B.  Gas  temperature 
variation  with  time  method 

Model  initial  temperature  is 
uniform. 

Required 

Preferred.  Flowever,  could  be 
analysed  if  the  variation  is 
known. 

Gas  temperature  distribution 
through  the  duct  must  be 
determined  at  the  start  of  the 
experiment,  before  the  surface 
temperature  rise  is  significant. 

Required 

No.  Flowever,  gas  temperature 
sensor  time  constants  must  be 
small  or  known  and  accounted 
for. 

The  relationship  between  the 
temperatures  measured  at  any 
plane  and  the  mixed  bulk 
temperature  ,(f),  is  known. 

Required 

Required 

Thermal  properties  of  substrate 
must  be  known. 

Not  required 

Yes 

Model  material  must  be 
homogenous 

No.  Only  needs  to 
have  uniform  initial 
temp. 

Yes.  Otherwise  analysis  is 
very  complicated. 

Table  1  Summary  of  the  requirements  for  the  2  methods. 


Gas  Temperature  Instrumentation 

Different  sensors  can  be  used  to  measure  the  gas  temperature.  These  include :- 

•  An  adiabatic  (nylon)  mesh  threaded  with  thermocouples.  This  approach  is  used  at 
Oxford-  see  Tsang  et  al.  (2000)  and  Mee  et  al.  (2000).  Figure  2  is  an  example  of  an 
arrangement  used  in  the  ribbed  duct  heat  transfer  research  discussed  below 

•  An  adiabatic  (nylon)  mesh  coated  with  liquid  crystals  then  monitored  with  a  video 
camera.  This  method  can  provide  a  detailed  temperature  profile  although  the 
disadvantage  is  that  optical  access  is  required.  See  Wang  et  al.  (1998)  and  Mee  et  al. 
(2000)  for  examples  of  its  application. 

•  An  adiabatic  (nylon)  mesh  threaded  with  platinum  resistance  thermometer(s).  This 
approach  produces  a  single  signal  for  the  temperature  averaged  along  the  length  of  the 
wire.  A  platinum  wire  threaded  through  the  mesh  has  a  resistance  change  proportional  to 
the  average  temperature  change  of  the  wire. 
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Figure  2  Thermocouples  woven  in  to  nylon 
mesh. 


Figure  3  Time  required  for  substrate 
dimensionless  temperature  to  change  5% 
for  different  materials  as  a  function  of  heat 
transfer  coefficient. 


Whichever  way  the  temperature  data  is  to  be  analysed,  the  experiment  needs  to  give  the  local 
mixed  bulk  (or  mixing  cup)  temperature  at  positions  within  the  model.  Evaluation  of  the 
mixed  bulk  temperature  at  any  plane  requires  an  integration  of  the  temperature  velocity 
product.  In  cases  where  it  is  not  possible  to  measure  the  velocity  distribution,  it  is  proposed 
that  the  latter  be  estimated  using  knowledge  of  the  total  coolant  flow  at  each  plane  and 
assumed  velocity  profiles.  Dependent  on  rig  size,  the  temperature  distribution  can  be 
determined  with  lower  uncertainty  using  several  temperature  sensors  at  each  plane.  The 
implication  of  uncertainty  in  the  profiles  on  uncertainty  in  h  is  discussed  in  detail 
below.  The  analysis  also  leads  to  guidance  on  the  minimum  separation  between  the  planes. 


The  link  between  the  measured  temperatures  and  the  local  mixed  bulk  temperature  is 
considered  by  introducing  a  profile  parameter,  (]),  that  relates  the  measured  gas  temperature, 
indicated  here  by  TExp,  to  the  true  mixed  bulk  temperature,  Tmb. 


4>i  = 


(3) 


The  subscript  i  refers  to  a  number  which  identifies  the  measurement  plane.  Note  that  the 
measured  temperature  can  be  the  average  temperature,  measured  using  for  example  by  a  PRT, 
or  can  be  can  be  determined  from  multiple  thermocouples.  Tw  is  the  model  wall  temperature. 


For  method  (A),  equation  (2)  requires  the  mixed  bulk  temperature  to  be  measured  and  so  the 
temperature  profile  factor  can  be  used  to  express  the  mixed  bulk  temperatures  in  terms  of  the 
measured  temperature.  Equation  (2)  becomes 

Qi- 1  ^exp(i-l) 

4>,  '  Texpi-Tw  \ 
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Evaluation  of  the  methods 

The  above  methods  were  used  to  determine  the  average  heat  transfer  coefficient  in  a  cooling 
passage  in  which  ribs  had  been  included  to  increase  the  heat  transfer.  The  research  was  part 
of  a  program  to  evaluate  the  performance  of  turbulators  with  a  larger  size  (relative  to  the 
passage  hydraulic  diameter)  than  those  traditionally  used  in  cooling  passages.  Some  of  the 
findings  of  the  work  are  presented  in  Tsang  et  al.  (2001).  Data  from  the  above  technique  are 
compared  below  to  those  evaluated  using  the  conventional,  transient  liquid  crystal  method. 
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Figure  4  Schematic  diagram  of  the  passage  used  for  rib  heat  transfer  research. 


Figure  4  shows  the  rig  together  with  the  positions  of  the  gas  temperature  measurement  planes. 
The  Perspex  rig  has  internal  section  with  50mm  square  sides.  The  air  is  drawn  from 
atmosphere  through  the  bell  mouth  intake  and  passes  immediately  through  a  fine  wire,  heater 
mesh  supplied  with  AC  power  from  a  transformer.  The  airflow  is  initiated  and  adjusted  to 
achieve  the  test  Reynolds  number  and  heat  transfer  starts  when  power  is  suddenly  switched  to 
the  mesh  heater.  Gas  temperature  signals  logged  throughout  a  typical  experiment  are  shown 
in  Figure  5. 


Figure  5  Logged  gas  temperatures  for  a  typical  experiment. 


At  each  plane,  the  gas  temperature  was  measured  using  fine  (76pm  wire)  diameter  Copper 
Constantan  thermocouples  threaded  through  an  adiabatic  mesh  as  described  in  Mee  et  al. 
(1999).  The  thermocouple  time  constant  reduces  as  the  flow  speed  increases.  The  former  was 
calculated  to  be  50ms  at  the  lowest  Reynolds  number  tested.  The  longest  flushing  time  was 
used  to  estimate  the  time  necessary  for  the  flow  field  to  establish.  Dividing  3m  by  the  lowest 
velocity  gives  a  flushing  time  of  less  than  1  second.  To  minimise  uncertainty  in  the  start 
temperature,  the  gas  temperatures  analysed  were  determined  by  curve  fitting  to  the  signals 
between  1  and  10  seconds  after  the  heat  initiation,  Figure  6.  Figure  2  shows  the  thermocouple 
positions  in  each  of  the  gas  temperature  measurement  planes.  All  of  the  thermocouples  were 
logged,  but  two  different  sets  of  thermocouples  were  used  in  the  analysis. 
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Figure  6  Extrapolation  used  to  determine  the  initial  gas  temperatures. 


Temperature  profile  factors 

The  profile  factor  was  evaluated  using  velocity  measurements,  Figure  7,  made  with  the 
miniature  probe  described  in  Tsang  et  al.  (2001).  Since  the  flow  through  the  smooth 
development  section  was  turbulent,  the  velocity  and  temperature  profiles  will  be  similar  by 
Reynolds  analogy.  This  means  that  the  distribution  of  the  dimensionless  velocity  (u/uc])  is  the 
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Figure  7  Velocity  distribution  at  plane  2.  This  is  downstream  of  the  smooth  development 
section  and  ahead  of  the  ribbed  walls. 


same  as  the  distribution  of  dimensionless  temperature  (T-Tw)/(T-Tcl). 
temperature  is  defined  as 


which  simplifies  to 


The  mixed  bulk 


(5) 


(6) 


since  the  density  and  specific  heat  are  uniform  across  the  measurement  plane.  A  little  algebra 
leads  to 


\ulud.(T-Tw)KTd-Tw)AA 
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and  Reynolds  analogy  gives 
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An  other  option  is  to  take  temperature  measurements  at  two  locations  which  have  the  same 
profile  factor.  The  two  equal  values  of  (|>  cancel  out  in  equation  (4).  An  example  would  be  to 
use  centreline  thermocouples  in  a  long  ribbed  passage  which  has  achieved  developed 
conditions.  The  temperature  profile  options  are  summarised  in  Table  2. 


Option 

Example 

Comments 

Profile  measured  directly  at  plane 

None 

Most 

accurate 

Profile  factor  inferred  from  velocity  measurement  and 
Reynolds  analogy. 

Rib  results 
present  papers 

Profile  shape  estimated  from  insight  into  the  flow  field. 

Eg.  the  flow  can  be  well  mixed  to  the  extent  that  the 
profile  can  safely  assumed  to  be  uniform. 

Impingement 
work  of  Son  et 
al.  (2001) 

Well  mixed 
profile  at 
exit  to  array. 

Analysis  applied  between  planes  for  which  the  profile  is 
identical. 

Developed 
flow  in  ribbed 
passage. 

Table  2  Summary  of  Profile  factor  evaluation  options. 


1.  Analysis  using  average  temperature  and  the  initial  temperature  method. 

The  above  integral  was  evaluated  at  plane  2  (Figure  4)  as  0.85  over  a  grid  of  2500 
measurement  points.  This  was  very  close  to  the  average  dimensionless  temperature 
difference  determined  at  the  9  thermocouple  measurement  sites.  The  latter  was  again 
determined  from  the  velocity  measurements  using  Reynolds  analogy  and  was  found  to  0.84. 
In  other  words,  the  average  gas  temperature  at  plane  2  was  an  excellent  approximation  to  the 
mixed  bulk  temperature. 


Figure  8  Ratio  of  heat  transfer  coefficient  determined  using  the  new  technique  (A)  to  h 
determined  from  a  conventional  liquid  crystal  heat  transfer  experiment. 

For  this  reason,  the  profile  factor  was  set  to  unity  at  the  inlet  and  the  average  temperature 
used  as  the  upstream  experimental  temperature  in  equation  (4).  All  of  the  rib  configurations 
tested  act  to  stir  the  flow  and  the  exit  temperature  profiles  were  found  to  be  sufficiently  flat 
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that  the  set  average  temperature  was  insensitive  to  the  set  of  thermocouples.  The  exit  gas 
temperature  profile  was  also  set  to  unity.  The  heat  transfer  coefficient  enhancement  values 
were  then  compared  to  data  determined  using  more  conventional  liquid  crystal  heat  transient 
experiments  and  the  ratio  between  the  two  values  compared  for  the  rib  tests.  The  results  are 
shown  in  Figure  8. 

2.  Analysis  using  core  representative  temperatures  and  the  initial  temperature 
method. 

In  some  situations,  it  will  not  be  possible  to  determine  an  average  temperature  representative 
of  the  mixed  bulk  temperature.  The  data  was  reprocessed  using  the  3  core  thermocouples  (C, 
F  and  G)  in  Figure  2  at  plane  2.  A  profile  factor  of  0.85  was  used  in  the  analysis  and  results 
(divided  by  the  normal  liquid  crystal  data)  are  presented  in  Figure  8. 


Uncertainty  analysis:  A.  Initial  temperature  distribution  method 

Uncertainty  analysis  shows  that  the  fractional  uncertainty  in  h  (assuming  that  uncertainties  in 
A,  m  and  cp  are  negligible),  is 


where  a  denotes  the  standard  deviation  in  a  parameter.  Recall  that  4>  is  a  factor  that 
characterises  the  temperature  profile.  As  defined  above,  it  is  the  quotient  of  (  mixed  bulk 
temperature  -wall  temperature)  divided  by  (measured  gas  temperature  minus  wall 
temperature).  The  ratio  cr^i/tft  is  the  fractional  uncertainty  in  this  temperature  profile 
parameter  at  the  plane  i.  The  equation  for  oh/h  shows  that  the  accuracy  in  the  measured  h  is 
reduced  if  the  gas  measurement  planes  are  too  close  together  since  the  NTU  between  planes 
reduces.  Figure  9  shows  the  calculated  uncertainty  for  different  uncertainty  levels  in  the 
temperature  profile  factors.  The  legend  indicates  the  permutation  of  inlet  and  outlet  profile 
factor  (  and  respective  uncertainties)  used.  Unsurprisingly,  the  accuracy  reduces  as  the  latter 
uncertainties  increase.  The  calculations  include  a  contribution  from  uncertainty  in  the 
measured  experimental  temperatures.  The  gas  temperatures  were  all  measured  with  a  data 
logger  fitted  to  a  PC  that  included  a  multiplexed  amplifier  and  cold  junction  compensation. 
Since  the  same  board  is  used  to  measure  the  gas  temperature,  the  uncertainty  in  the  gas 
temperature  differences  from  the  stalling  temperature  is  considerably  lower  than  uncertainty 
in  the  absolute  temperature.  A  level  of  0.3°C  for  the  standard  error  in  temperature  difference 
for  both  planes  was  used  in  Figure  9.  The  exit  plane  temperature  difference  is  fixed  by  the 
inlet  plane  temperature  difference  and  the  NTU  between  planes  so  the  above  equation 
becomes 


1 


h  NTU i 


'rexp-v^_, 


^  -rJexpt-ATT/) 
Vi 


(10) 


From  the  chart  it  is  clear  that  the  measurement  planes  need  to  be  sufficiently  far  apart  to 
achieve  an  acceptable  accuracy.  The  uncertainty  stalls  to  rise  as  the  NTU  goes  beyond  2 
(dependent  on  the  parameters)  since  the  last  term  in  equation  (10)  stalls  to  dominate. 
Physically,  the  NTU  needs  to  be  sufficiently  large  between  planes  to  ensure  the  gas 
temperature  difference  is  measurable,  but  not  so  large  that  the  exit  plane  temperature  is  the 
same  as  the  wall  temperature. 

The  NTU  term  in  the  denominator  of  equation  (10)  is  responsible  for  the  initial  reduction  of 
uncertainty  with  NTU.  It  is  pleasing  to  note  that  the  uncertainty  in  h  can  be  less  than  the 
uncertainty  in  the  profile  factors. 
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Figure  9  Calculated  uncertainty  in  h  as  a  function  of  the  Number  of  Transfer  Units  between 
planes  for  different  shape  factor  uncertainty  levels.  The  uncertainty  in  both  gas  temperature 
minus  initial  temperature  is  set  at  0.3°C  and  the  inlet  gas  to  starting  temperature  difference  is 
40°C. 
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Figure  10  The  diamond  symbols  show  the  calculated  uncertainty  in  h  as  a  function  of  the 
Number  of  Transfer  Units  between  planes  for  method  B.  The  uncertainty  in  both  gas 
temperature  minus  initial  temperature  is  set  at  0.3°C  and  the  uncertainty  in  the  profile  factors 
is  10%.  The  legend  indicates  the  initial  difference  between  the  upstream  gas  temp  and  the 
model. 

Uncertainty  analysis:  B.  Gas  temperature  variation  with  time  method 

A  Monte-Carlo  technique  was  used  to  determine  the  uncertainty  in  the  h  measurement  for  the 
second  method  and  the  results  for  different  NTUs  are  plotted  in  Figure  10.  The  uncertainty 
shows  the  same  trends  as  before  and  has  approximately  the  same  level. 
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Conclusions 

A  new  way  of  processing  gas  temperature  data  from  apparatus  used  for  transient  liquid  crystal 
experiment  has  been  described  and  evaluated.  The  data  resolution  is  reduced,  relative  to  a 
liquid  crystal  experiment,  but  the  method  offers  advantages  in  situations  where  processing 
time  is  restricted  or  optical  access  is  not  possible.  The  experimental  uncertainty  has  been 
evaluated  and  acceptable  accuracy  is  achieved  when  the  there  is  sufficient  temperature  change 
between  measurement  planes.  The  new  method  offers  the  potential  to  perform  internal 
cooling  experiments  in  metallic  models  at  with  higher  heat  flux. 
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Nomenclature 

A  area,  nr 

cp  specific  heat  at  constant  pressure,  Jkg  ’K1 

h  heat  transfer  coefficient,  Wrn  2  K1 

m  mass  flow  rate,  kgs  1 

NTU  number  of  transfer  units 

q  heat  flux,  Wm 2 

T  temperature,  °C 

Subscripts 

exp  measured  temperature 

exp-W  difference  between  measured  and  stalling  wall  temperatures 
i  plane  number 

w  initial  wall 

Greek 

cf)  temperature  profile  factor 

p  density,  kgm 3 
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Paper  Number:  4 

Name  of  Discusser:  J.  von  Wolfersdorf,  ABB  Baden  Switzerland 

Question: 

How  did  you  measure  the  heat  transfer  at  the  high  blockage  ribs  in  your  conventional  liquid 
crystal  technique? 

Answer: 

The  ribs  were  mode  from  copper  and  were  fixed  atop  of  the  liquid  crystal  coating.  The 
conductivity  of  the  rib  is  sufficiently  high  that  the  analysis  for  a  calori  meter  ontop  of  a  semi¬ 
infinite  substrate  can  be  used  to  determine  each  rib  heat  transfer  coefficient. 

The  method  is  described  in  Wang  et  al.  1996.  Heat  Transfer  measurements  to  a  gas  turbine 
cooling  passage  with  inclined  ribs.  Journal  of  Turbomachinery,  Vol.  170.  pp  63-69. 
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